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ABSTRACT 

Context. We report on the results of new simulations of near-infrared (NIR) observations of the Sagittarius A* (Sgr A*) 
counterpart associated with the super-massive black hole at the Galactic Center. 

Aims. Our goal is to investigate and understand the physical processes behind the variability associated with the NIR 
flaring emission from Sgr A*. 

Methods. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European 
Southern Observatory's Very Large Telescope and CIAO NIR camera on the Subaru telescope (13 June 2004, 30 July 
2005, 1 June 2006, 15 May 2007, 17 May 2007 and 28 May 2008). We used a model of synchrotron emission from 
relativistic electrons in the inner parts of an accretion disk. The relativistic simulations have been carried out using the 
Karas-Yaqoob (KY) ray-tracing code. 

Results. We probe the existence of a correlation between the modulations of the observed flux density light curves 
and changes in polarimetric data. Furthermore, we confirm that the same correlation is also predicted by the hot spot 
model. Correlations between intensity and polarimetric parameters of the observed light curves as well as a comparison 
of predicted and observed light curve features through a pattern recognition algorithm result in the detection of a 
signature of orbiting matter under the influence of strong gravity. This pattern is detected statistically significant 
against randomly polarized red noise. Expected results from future observations of VLT interferometry like GRAVITY 
experiment are also discussed. 

Conclusions. The observed correlations between flux modulations and changes in linear polarization degree and angle 
can be a sign that the NIR flares have properties that are not expected from purely random red-noise. We find that 
the geometric shape of the emission region plays a major role in the predictions of the model. From fully relativistic 
simulations of a spiral shape emitting region, we conclude that the observed swings of the polarization angle during NIR 
fiares support the idea of compact orbiting spots instead of extended patterns. The effects of gravitational shearing, 
fast synchrotron cooling of the components and confusion from a variable accretion disk have been taken into account. 
Simulated centroids of NIR images lead us to the conclusion that a clear observation of the position wander of the 
center of NIR images with future infrared interferometers will prove the existence of orbiting hot spots in the vicinity 
of our Galactic super-massive black hole. 

Key words, black hole physics: general, infrared: general, accretion, accretion disks. Galaxy: center. Galaxy: nucleus 



1. Introduction 

The nearest super-massive black hole candidate 4 x 
10^ Mq) lies at the center of our galaxy, as inferred 
from motions of stars near the Galactic Center (Eckart 
& Genzel 1996, 1997; Eckart et al. 2002; Schodel et al. 
2002; Eisenhauer et al. 2003; Ghez et al. 2000, 2005, 
2008; Gillessen et al. 2009). With a luminosity of 10"^ - 
10~^^ LEdd^ where LEdd is its limiting Eddington lumi- 
nosity, Sagittarius A*, the radio source associated with 
this SMBH, is one of the most extreme sub-Eddington 
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sources accessible to observations. However, X-ray and 
near-infrared (NIR) flares are routinely detected with high 
spatial and spectral resolution observations (Baganoff et al. 
2001; Porquet et al. 2003, 2008; Genzel et al. 2003; Eckart 
et al. 2004, 2006a-c, 2008a-c; Meyer et al. 2006a,b, 2007; 
Yusef-Zadeh et al. 2006a,b, 2007, 2008). These short bursts 
of increased radiation last normally for about 100 minutes 
and occur four to five times a day (see Fig. 1 for a typical 
behavior of Sgr A* during a flaring phase in NIR band- 
width). 

Recent NIR and X-ray observations have revealed the 
non-thermal nature of high frequency radiation from Sgr A* 
(Eckart et al. 2006a-c, 2008a-c; Gillessen et al. 2006; 
Hornstein et al. 2007). Sgr A* is probably visible in the 
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Fig. 1. Sgr A* as it was observed in NIR L'-band (3.8 fim) on 3 June 2008 between 05:29:00 - 09:42:00 (UT time), {a-f) 
show the observed images of Sgr A* after 0, 50, 91, 135, 155 and 212 minutes off the start of observation. The images 
show that when Sgr A* is in its flaring state the flux changes up to 100% in time intervals of the order of only tens of 
minutes. 



NIR regime only during its flaring state. The short time 
scale variabilities seen during several observed NIR and 
X-ray flares argue for an emitting region not bigger than 
about ten Schwarzschild radii (r^ = = 2r^ = 1.2 x 

( 4x106 Mq ) associated super-massive black 

hole (Baganoff et al. 2001; Genzel et al. 2003). We have 
scaled the relevant physical distances according to the grav- 
itational radius (vg) throughout this paper. The NIR flares 
are highly polarized and normally have X-ray counterparts, 
which strongly suggests a synchrotron-self-Compton (SSC) 
or inverse Compton emission as the responsible radiation 
mechanism (Eckart et al. 2004, 2006a,b; Yuan et al. 2004; 
Liu et al. 2006). Several observations have already con- 
firmed the existence of a time lag between the simultaneous 
NIR/X-ray flares and the flares in the lower frequencies. 
This is interpreted as a sign for cooling down via adiabatic 
expansion (Eckart et al. 2006a, 2008b, c; Yusef-Zadeh et al. 
2006a,b, 2007, 2008; Marrone et al. 2008, Zamaninasab et 
al. 2008a). 

The other feature related to these NIR/X-ray flares are 
the claimed quasi-periodic oscillations (QPOs) with a pe- 
riod of 20 ± 5 minutes, which have been reported in several 
of these events (Genzel et al. 2003; Belanger et al. 2006; 
Eckart et al. 2006b,c; 2008a; Meyer et al. 2006a,b, Hamaus 
et al. (2009)). Short periods of increased radiation (the so 
called "NIR flares", normally around 100 minutes) seem to 
be accompanied by QPOs. All the studies mentioned above 



probed this 20±5 minutes quasi-periodicity, by performing 
a sliding window analysis with window lengths of the order 
of the flaring time. Recently, Do et al. (2009) argued that 
they did not find any significant periodicity at any time 
scale while probing their sample of observations for a pe- 
riodic signal. Their method is based on the Lomb-Scargle 
periodogram analysis of a sample of six light curves and 
comparing them with several thousands of artificial light 
curves with the same underlying red-noise. One must note 
that the suggested QPOs are transient phenomena, last- 
ing for only very few cycles (50-100 minutes). This kind of 
behavior, along with the inevitable uncertainty in the red 
noise power law index determination, makes a clear and 
unambiguous detection (> 5cr) of a periodic signal very 
difficult. Whenever a flare of Sgr A* was observed with 
polarimetry, it was found that it is accompanied by signifi- 
cant polarization which varies on similarly short timescales 
as the light curve itself. By carrying out an analysis only on 
the total flux, some pieces of the observed information are 
ignored. It is already known that polarimetric data have 
been shown to be able to reveal substructure in flares, even 
when the light curve appears largely featureless (e.g. see 
Fig. 4 in Eckart et al. 2006b). The other main advantage 
of polarimetric observations is that, in addition to the fiux 
density light curve, one can analyze the changes in the ob- 
served degree of polarization and the changes in polariza- 
tion angle during flaring time as well. 
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Fig. 2. Our sample of light curves of Sgr A* flares observed in NIR Ks band (2.2/im) polarimetry mode. The events were 
observed on 2004 June 13 (a), 2005 July 30 (b), 2006 June 1 (c), 2007 May 15 (d), 2007 May 17 (e) and 2008 May 28 (f). 
In each panel, the top shows the de-reddened flux density measured in mJy (black), the middle shows the polarization 
angle (blue) and the bottom shows the degree of linear polarization (red) . The gaps in the light curves are due to the 
sky background measurements. 



The claimed quasi-periodicity has been interpreted as 
being related to the orbital time scale of the matter in 
the inner parts of the accretion disk. According to well- 
known observed high frequency quasi-periodic oscillations 
(HFQPOs) in X-ray light curves of stellar mass black holes 
and binaries (Nowak & Lehr 1998), this interpretation is of 
special interest since it shows a way to better understand 
the behavior of accretion disks for a wide range of black 
hole masses. The recent unambiguous discovery of a one 
hour (quasi-)periodicity in the X-ray emission light curve of 
the active galaxy RE J 1034+396 provides further support 
to this idea and extends the similarity between stellar-mass 
and super-massive black holes to a new territory (Gierlihski 
et al. 2008). 



Although the origin of the observed QPOs in the sources 
associated with black holes is still a matter of debate, sev- 
eral magnetohydrodynamic (MHD) simulations conflrmed 
that it could be related to instabilities in the inner parts of 
the accretion disks, very close to the marginally stable or- 
bit of the black hole (rmso), and also possibly connected 
with the so-called "stress edge" (Hawley 1991; Chan et 
al. 2009b). If the flux modulations are related to a sin- 
gle azimuthal compact over-dense region (hereafter: "hot 
spot"), orbiting with the same speed as the underlying 
accretion disk, one can constrain the spin of the black 
hole by connecting the observed time scales of QPOs to 
the orbital time scale of matter around the black hole: 
T = 2.07(ri + ^)( 4xi^m. ) ™^ (Bardeen et al. 1972, 
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Date (Telescope) 



Spectral domain UT start UT stop Max. flux Min. flux Average flux Average polarization 
time time sampling rate sampling rate 



2004 June 13 (NACO) 


2.2/im 


07:20:02 


09:15:08 


5.19 mJy 


2.21 mJy 


1.2 min 


3 min 


2005 July 30 (NACO) 


2.2/im 


02:07:36 


07:03:39 


8.19 mJy 


1.23 mJy 


1.2 min 


3 min 


2006 June 1 (NACO) 


2.2/im 


04:26:03 


10:44:27 


19.33 mJy 


0.72 mJy 


1.5 min 


2 min 


2007 May 15 (NACO) 


2.2/im 


09:08:14 


09:42:12 


22.27 mJy 


1.70 mJy 


1.5 min 


2 min 


2007 May 17 (NACO) 


2.2/im 


04:42:14 


09:34:40 


11.87 mJy 


1.86 mJy 


1.5 min 


2 min 


2008 May 28 (CIAO) 


2.15/im 


09:22:51 


13:00:37 


7.70 mJy 


0.97 mJy 


3.3 min 


10 min 



Table 1. Observations log. 



where — 1 < a < 1 is the black hole dimensionless spin pa- 
rameter and r is the distance of the spot from the black 
hole). The characteristic behavior of general relativistic 
flux modulations produced via the orbiting hot spots have 
been discussed in several papers (see e.g. Cunningham & 
Bardeen 1973; Abramowicz et al. 1991; Karas & Bao 1992; 
Hollywood et al. 1995; Dovciak 2004, 2007). In this paper, 
we have used a spotted accretion disk scenario to model the 
observed patterns of our sample of NIR light curves. 

Several authors have proposed different models in or- 
der to explain the flaring activity of Sgr A*. These models 
cover a wide range of hypotheses like disk-star interactions 
(Nayakshin et al. 2004), stochastic acceleration of electrons 
in the inner region of the disk (Liu et al. 2006), sudden 
changes in the accretion rate of the black hole (Liu et al. 
2002), heating of electrons close to the core of a jet (Markoff 
et al. 2001; Yuan et al. 2002), trapped oscillatory modes in 
the inner regions of the accretion disk in the form of spiral 
patterns or Rossby waves (Tagger & Melia 2006; Falanga et 
al. 2007; Karas et al. 2008), non-axisymmetric density per- 
turbations which emerge as the disk evolves in time (Chan 
et al. 2009b), non-Keplerian orbiting spots falling inward in- 
side the plunging region created via magnetic reconnections 
(Falanga et al. 2008), and also comet like objects trapped 
and tidally disrupted by the black hole (Cadez et al. 2006; 
Kostic et al. 2009). 

Observational data render some of these models un- 
likely. The star-disk interaction model is unable to produce 
the repeated flux modulations and also the high degree of 
polarization since it mainly deals with thermal emission. 
Tidal disruption of comet-like objects are also unable to 
reproduce the observed rate of flares per day, since the es- 
timated capture rate of such objects for the Sgr A* environ- 
ment is at least one order of magnitude lower (Cadez et al. 
2006). Nevertheless, several viable models exist and make 
different predictions that can be distinguished observation- 
ally. For example, one important characteristic prediction 
of hot spot models is about the wobbling of the center of the 
images (Broderick & Loeb 2006a, b; Paumard et al. 2006; 
Zamaninasab et al. 2008b; Hamaus et al. 2009). Significant 



effort has been already devoted to measure this possible 
position wander of Sgr A* in the mm, sub-mm and NIR 
regimes (Eisenhauer 2005b, 2008, Gillesen 2006, Reid et al. 
2008, Doelleman et al. 2008). In this paper we discuss how 
NIR polarimetry and the next generations of VLT inter- 
ferometry (VLTI) and Very Long Baseline Interferometry 
(VLBI) experiments can provide data to support or re- 
ject certain models for the accretion flow/outflow related 
to Sgr A*. Obtaining accurate data on the accretion flow 
of Sgr A* can lead us to a better understanding of the 
physics of strong gravitational regimes, formation of black 
holes and their possible relation to the galaxy formation 
process in a cosmological context. 

In Sect. 2 we present a complete sample of NIR light 
curves observed in the polarimetry mode. A brief descrip- 
tion about the details of the observation and data reduc- 
tion methods is provided. We discuss the quasi-periodicity 
detection methods and present the results of a correlation 
analysis between the flux and polarimetric parameters. A 
general description of our model setup and results of simu- 
lations are discussed in Sect. 3. We show how NIR polarime- 
try can be used as a way to constrain physical parameters 
of the emitting region (like its geometrical shape) in Sect. 
4 and Sect. 5. In Sect. 6 we mainly discuss the predictions 
that the future NIR interferometer (GRAVITY) is expected 
to reveal and how different assumptions in the model pa- 
rameters can modify the results. In Sect. 7 we summarize 
the main results of the paper and draw our conclusions. 

2. NIR polarimetry 

2.1. Observations and data reduction 

All observations we refer to in this paper have been carried 
out in the NIR Ks band with the NIR camera CONICA 
and the adaptive optics (AO) module NAOS (NACO) at 
the ESO VLT unit telescope 4 (YEPUN) on Paranal, Chil^ 

^ Based on observations at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO) on Paranal in Chile; 
Programs:075.B-0093 and 27LB-5019(A). 
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and CI AO NIR camera on the Subaru telescope These fa- 
cilities are suited for both time resolved observations of to- 
tal intensity and polarimetric degree and angle with a sam- 
pling of about two to three minutes. The N AOS/CONIC A 
NIR camera installed on UT4, VLT allows for a simulta- 
neous measurement of two orthogonal directions of electric 
field vector via a Wollaston prism. The combination with a 
rotary half-wave plate allows the rapid alternation between 
measurements of different angles of the electric vector. This 
is crucial for determining the linear polarization character- 
istics of a time- varying source. The CI AO camera uses a 
rotating half-wave plane combined with a fixed wire grid 
polarizer for measuring the linear polarization. 

Since the first NIR polarimetric observation of Sgr A* 
in 2004 (Eckart et al. 2006a), several polarized flares have 
been observed (Meyer et al. 2006 a,b, 2007; Eckart et al. 
2008a). In the VLT observations, the infrared wavefront 
sensor of NAOS was used to lock the AO loop on the 
NIR bright (K-band magnitude ~ 6.5) supergiant IRS 7, 
located about 5.6^' north of Sgr A*. Atmospheric condi- 
tions (and consequently the AO correction) were stable 
enough during the observations for doing high angular res- 
olution photometry and polarimetry (with a typical coher- 
ence time of 2 milliseconds and larger). The exposures have 
been dithered. The reductions of the 1 June 2006 and 15 
May 2007 data presented here have been repeated for this 
publication to confirm the significance of the discussed fea- 
tures (see also Eckart et al. 2008a). The 17 May 2007 data 
have not been published before. All exposures were sky sub- 
tracted, flat-fielded, and corrected for dead or bad pixels. 
As the most important improvement of the new reduction 
the dithered exposures have been aligned with sub-pixel 
accuracy by a cross-correlation method (Devillard 1999). 
PSFs were extracted from these images with StarFinder 
(Diolaiti et al. 2000). The images were deconvolved with 
the Lucy- Richardson (LR) algorithm. Beam restoration was 
carried out with a Gaussian beam of FWHM corresponding 
to the respective wavelength. The final resolution at 2.2/im 
is about 60 milli-arcseconds (mas). Flux densities of the 
sources were measured by aperture photometry. Because of 
the high accuracy of the image alignment it was possible 
to separate Sgr A* from the nearby stars S17 and S13 by 
choosing a circular aperture of about 52 mas radius (see 
Fig. 15 of Eckart et al. 2006a and discussion therein), re- 
sulting in a better correction for the flux contribution of 
these stars. The data was corrected for extinction using 
= 2.8 (Eisenhauer et al. 2005a, Schodel et al. 2007). 
The relative flux density calibration was carried out using 
known K-band flux densities and positions of 14 sources in 
the IRS16 cluster by R. Schodel (private communications). 
This results in a K-band flux of the high velocity star S2 of 
22 zb 1 mJy, which compares well with the magnitudes and 
fluxes for S2 quoted by Chez et al. (2005b) and Genzel et 
al. (2003). The measurement uncertainties for Sgr A* were 
obtained from the reference star S2. For more details about 
CIAO observations see Nishiyama et al. (2009). 



^ Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 



2.2. Data analysis 
2.2.1. Periodicity 

Figure [2] shows flare events observed in the NIR K-band 
(2.2/im) on 13 June 2004, 30 July 2005, 1 June 2006, 15 
May 2007, 17 May 2007 and 28 May 2008. This sample in- 
cludes all flare events observed in NIR polarimetry during 
recent years according to the knowledge of the authors. The 
flux densities rise and come back to their quiescent level in 
time intervals of roughly 100 minutes. The measured val- 
ues of total flux, degree of polarization, and angle of po- 
larization vary significantly on time scales of ~10 minutes. 
These abrupt changes can be more clearly detected when 
the events are in their brightest state. Even if the flares are 
different in some aspects (e.g. the ratio of the changes of 
the flux, the maximum brightness achieved, or the degree 
of linear polarization), there could exist some features that 
are repeated in our sample. Here we perform a quantitative 
analysis in order to probe such features. We focus first on 
detecting periodic signatures in flux densities and then a 
possible correlation between changes in flux and polarimet- 
ric light curves. 

The autocorrelation function and Lomb-Scargle peri- 
odograms can be used to detect signatures of time peri- 
odic structures. For the autocorrelation analysis we used 
the z-transformed discrete correlation function (ZDCF) al- 
gorithm (Alexander 1997), which is particularly useful for 
analyzing sparse, unevenly sampled light curves. For Lomb- 
Scargle periodograms we followed Press and Rybicki (1989). 

Figure [3] shows the cross-correlation of the flux density 
light curves of our sample with themselves by using ZDCF 
method. The ZDCF of the 30 July 2005, 1 June 2006, 15 
and 17 May 2007 flares show peaks around 20 ±5 min time- 
lag, and specially for 15 and 17 May 2007 the peaks look 
significant. These peaks can be signs of a possible period- 
icity. 

As one can see, only 17 May 2007's ZDCF shows a peak 
above the 99.8% significance threshold. The false alarm val- 
ues have been derived by repeating the same ZDCF anal- 
ysis on the 10^ random red-noise light curves. The criteria 
for this comparison and the way the light curves are pro- 
duced are described below. Here we must note that since 
the ZDCF values are normalized and distributed in a [-1:1] 
interval, the distribution of the values for each time lag (r) 
is not Gaussian (see Fig. |4|. As a result we have used the 
median and percentile nomenculture instead of the normal 
standard deviation formalism for deriving the false alarm 
levels. 

We have highlighted the time windows in which ZDCFs 
show a peak around 20 ±5 min, coinciding with the claimed 
quasi-periodicity. Even though the 1 June 2006 and 15 May 
2007's ZDCFs do not show peaks above 99.8% significance, 
they however still reach 99% and 95.4% levels, respectively. 
Also, by a rough estimate, the significance of the ZDCF 
peaks are correlated to the brightness of the flares, and for 
the faint events, as is observed for example in 2004, there 
is no detectable peak. 

Figure |5] shows the Lomb-Scargle periodogram for each 
flux light curve. The usual factor of four times over- 
sampling was applied in order to increase the sensitivity. All 
periodograms show the power law behavior, V (x , show- 
ing greater power at lower frequencies. Here V is power, / is 
frequency and u is power law index. All light curves follow 
a red noise under-lying power spectral density (PSD) with 
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Fig. 3. ZDCF of the flux light curves of 13 June 2004 (a), 30 July 2005 (b), 1 June 2006 (c), 15 May 2007 (d), 17 May 
2007 (e) and 28 May 2008 (f). The vertical colored boxes indicate the position of the closest peaks to the zero time lags. 
Dotted lines show the median (red), 68.3% (green) and 99.8% (blue) thresholds derived from 10^ red noise simulated 
light curves. 




the periodogram we have repeated the Lomb-Scargle anal- 
ysis for 10^ simulated red noise light curves. For simulating 
the red noise we have followed the algorithm by Trimmer 



Fig. 4. ZDCF of 10^ simulated red noise light curves over- 
plotted in one image (a). The resulted average and standard 
deviation as well as the median and percentile values are 
presented. The distribution of the ZDCF values for differ- 
ent time lags (r) are also depicted (b-e) , which clearly show 
non-Gaussian distributions, specially for small values of r. 



-1.19 < CO < —1.77. To test the significance of the peaks in 



IT] for examples). The 
or each event sepa- 



and Konig (1995) (see Figs. [9) [iq and 
red noise light curves were producec. 
rately following the procedure below: (1) Determining the 
slope of the observed PSD for each event using a linear fit 
to the periodogram in log-log space and also by fitting a 
first-order autoregressive function to the PSD (Schulz and 
Mudelsee 2002) and averaging the results of both methods 
(the resulted co is presented in the lower right corner of each 
plot in Fig.|5|. (2) 10"^ red noise light curves produced with 
the related co following the method of Timmer and Konig 
(1995). The light curves have been produced by selecting 
a middle part of a light curve with a length at least ten 
times longer than the observed light curve following the ar- 
gument by Uttley et al. (2002). The selected segment of 
the light curve is sampled to the same time bin of the cor- 
responding observation. This will correct for any artificial 
effects caused by time lags which exist on our light curves 
and also the uneven sampling of our observations. (3) The 
value of the flux has multiplied by one factor as the mean of 
the simulated flux has the same value as the observational 
one for each night. 

None of the periodograms of Fig. [5] show a peak that ex- 
ceeds the 99.8% threshold except the ones correlated to the 
lengths of the flare events (~100 minutes). We specially 
highlighted the windows in which corresponding ZDCFs 
peak. In order to see if there exists a persistent frequency 
peak in the PSDs, we averaged all observed periodograms 
(Fig. |6] (b)). The average periodogram doesn't show a sig- 
nificant peak in comparison with the false alarm level even 
if we exclude the 13 June 2004 light curve (Fig. |6] (c)). 
As we mentioned before the short time scale of the flares 
and limited number of cycles make a significant detection 
of any periodicity very difficult. In addition, the possible 
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Fig. 5. Lomb-Scargle periodograms of the flux light curves for 13 June 2004 (a), 30 July 2005 (b), 1 June 2006 (c), 15 
May 2007 (d), 17 May 2007 (e) and 28 May 2008 (f). The dashed lines show the median (red) and , 68.3% (green) and 
99.8% (blue) thresholds derived from the 10^ red noise simulated light curves. The highlighted boxes show the band of 
frequencies that have corresponding ZDCF peaks. 




[l/mln] 



Fig. 6. Lomb-Scargle periodograms of the flux light curves 
of our sample overplotted in the same plot (a). Averaged 
result of all periodograms (b) and excluding the 13 June 
2004 flare (c). Highlighted box shows the expected 20 ± 5 
minutes periodicity band. 



quasi-periodic structure can vary during the time of the 
flare if it is connected to the falling clumps of matter into 
the black hole (Falanga et al. 2008). In this case the Lomb- 




f [l/min] 



Fig. 7. Same as Fig. [e] but averaged for a band of frequen- 
cies with the same size as the highlighted region. 



Scargle algorithm finds different frequencies and allocates 
them separate values of power. This effect can result in a 
periodogram in which the values of power are higher in a 
band of frequencies instead of a specific value. As one can 
see in Fig.[5j 1 June 2006 and 15 May 2007 can be the candi- 
dates for such an effect since the value of the PSD function 
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Fig. 8. Cross-correlation between the flux and polarization angle (degree) light curves of Fig.|2|[13 June 2004 (a), 30 
July 2005 (b), 1 June 2006 (c), 15 May 2007 (d), 17 May 2007 (e) and 28 May 2008 (f)]. In each panel, top (bottom) 
shows the correlation between the fiux and the polarization angle (degree of linear polarization). Dashed lines indicate 
the position of the median (red), 68.3% (green) and 99.8% (blue) significance levels derived from 10^ simulated red noise 
light curves. 



remains high in the highlighted window. To test for such 
an effect we averaged the periodograms with a new bin size 
and repeated the same procedure for the random red noise 
PSDs (Fig. [7|). The resulting averaged periodogram does 
not show any significant peak again (Fig. It] (b)), even if 
we exclude the 13 June 2004 event (Fig. m (c)). Since this 
procedure is very sensitive to the size oithe chosen win- 
dow (which can vary from event to event), a more detailed 
analysis is needed to study the possible evolution of any 
periodic signal. 



2.2.2. Importance of polarimetry 

We have also performed cross correlation analysis between 
variations of flux and degree (angle) of polarization. A 
search for any short time-lag correlation has been carried 
out by scanning the light curves using a sliding window 
method (see Fig. [12] panel (a)). The size of the scanning 
window (e) is fixecTon 40 minutes since we are interested in 
magnifying any short lag correlation related to the possible 
20 ± 5 minutes quasi-periodicity. Figure [S] shows the cross- 
correlations of the two sets of mentioned light curves for all 
events in our sample. The time steps of the scans were fixed 
to be five minutes. The presented results are the average of 
all scans weighted by the flare signal to noise ratio (i.e. each 
part of the flare that is brighter, has more weight). Since the 
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Fig. 9. (a): Simulated light curves of four different polari- 
metric channels all showing white noise behavior {cj = 0). 
(b): Total flux, (c): Angle of polarization, (d): Degree of 
polarization, (e): Lomb-Scargle periodograms of the total 
flux. 



ZDCF algorithm needs at least 11 points per bin for reliable 
results (specially for the error estimation; see Alexander 
1997) a linear interpolation of the polarimetric data points 
has been performed. We must mention that due to the con- 
fusion from nearby stars and the diffuse background emis- 
sion the accuracy of measuring polarimetric parameters is 
related to the brightness of Sgr A*. This means that the 
most reliable polarimetric data are measured when strong 
flares happen, specially when the source is in its brightest 
state. 

In order to test whether a random red noise model can 
produce the same correlation patterns we need to simulate 
red noise light curves including polarimetric data. For this 
purpose, we have simulated random E vectors for four per- 
pendicular directions: 




Time [min] 




0.10 



Eg = Eee 



(1) 



Fig. 10. Same as Fig. [9] for a; = — 1 (flicker noise). 



where ^ = 0°, 45°, 90°, 135°, is the corresponding normal 
vector, (j) is the phase and Eq is a random value following 
the algorithm by Trimmer and Konig (1995). Following the 
detection of the signal including the cross-talk from non- 
orthogonal neighboring channels, the polarized flux in each 
channel is then produced by using the Mueller matrix for- 
malism. The resultant four different channels (^0,^45,^90 
and F135) all show the same power law index (Figs. |9j [To] 
andpT] (a)). Using 

(2) 
(3) 
(4) 

(5) 
(6) 

(where F is the total flux, Q and U are the Stokes param- 
eters, X is the polarization angle and I!/, is the degree of 
linear polarization) the flux and polarization light curves 



F = 


^0 + ^90 


Q = 


^0 — ^90 


u = 


^45 — ^135 


X = 


-arctan(-) 


Hl = 




F 
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Fig. 11. Same as Fig. [ojfor cj = —2 (red noise). 



for each set can be derived. Figures [21 [To] and [11] show ex- 
amples of such light curves with three different PSD slopes: 
CO = (white noise), co = —1 (flicker noise) and co = —2 (red 
noise). 

As one can see in Fig. 11 for a; = —2 some random 
correlations between the changes in the total flux and the 
polarimetric parameters can occur. In order to see whether 
the observed correlations in our sample are just the same 
random coincidences or if they show signs of a more sub- 
tle process, we have simulated 10^ artificial light curves for 
each value of uj as derived from the average PSD of obser- 
vations. Then for each set of light curves the correlation 
function between flux and polarization angle (degree) has 
been calculated. Dashed lines in Fig. [8] show the median, 
68.3% and 99.8% false alarm values. One can see that some 
of the observed correlations are above 99.8% significance 
level. Even though most of the mentioned cross-correlation 
peaks happen around the same value (around zero time- 
lag) not all of them are exactly in the same r. To examine 
how probable it is that a strong deviation from the average 
in the red noise simulation repeatedly happens in a spe- 
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Fig. 12. (a) Sketch showing how the cross-correlations of 
Fig. [8] have been derived. A movins; window of the size e 
scans the flux and polarization angle (degree) light curves. 
The final cross-correlation is the average of all windows, 
(b) Sketch showing how the probability function in Fig. 
[13] has been derived. The algorithm finds the number of 
events in which two ZDCFs show at least one point over 3<j 
threshold of red noise correlation in a window the size of S. 
By repeating the same procedure for 10^ cross-correlation 
of red noise light curves the probability that a significant 
correlation is expected to happen in that window is derived. 



cific window we have calculated the probability that the 
mentioned cross correlations show significant peaks (above 
99.8%) in a window of a size 6 (see Fig. 12 (b)). 



Figure [13] shows this probability as derived for the simu- 
lated light curves (solid line) and observations (circles and 
triangles). This analysis shows that it is very unlikely to 
observe a correlation between total flux and polarimetric 
parameters approximately at the same time lags (small ^), 
while our observations show that a strong correlations ex- 
ist in the light curves of Sgr A* and that they also repeat 
themselves for approximately the same time lag. We must 
note here that the value derived for the observed proba- 
bility in Fig. [13] (circles and triangles) are derived from a 
sample of only six sets of light curves. In order to make a 
more reliable statistical analysis, more NIR observations of 
Sgr A* in polarimetric mode needs to be done in the fu- 
ture. Furthermore, the method described here can be used 
in principle for polarimetric observations of other sources 
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Fig. 13. Top: Probability that two sets of 10^ simulated 
red noise light curves show significant correlation between 
total flux and polarization angle in a window of the size of S 
minutes (solid line). Bottom: Same as top panel for the cor- 
relation between flux and polarization degree. Dotted lines 
show la (red), 3a (green) and 5<j (violet) levels of confi- 
dence for the probability function, calculated by repeating 
the whole analysis 100 times. The circles and triangles indi- 
cate the probabilities derived from our sample of NIR light 
curves (see Fig. [2]and|8|. 

showing the same variability; which may help in under- 
standing the general underlying physical process causing 
this kind of behavior. 

Without GR effects being taken into account, the phys- 
ical models which have been already proposed to simulate 
the observed red-noise light curves of AGNs (Lyubarskii 
1997; Armitage & Reynolds 2003, Vaughan et al. 2003) 
would have difficulties in reproducing this type of correla- 
tion between the behavior of polarimetric parameters and 
the total flux. As a result, the observed correlation be- 
tween changes in flux and the polarimetric data suggests 
a way to distinguish between the possible physical pro- 
cesses responsible for the overall red- noise behavior. A semi- 
analytical study by Pechacek et al. (2008) showed that a 
signal generated by an ensemble of spots randomly cre- 
ated on the accretion disk surface can produce red noise 
signals with PSD slopes of the order of -2. In their simula- 
tions the spot generation is governed by Poisson or Hawkes 
processes. In combination with our observations of Sgr A* 
the spotted disk scenario is a possible explanation for this 
commonly observed red noise behavior, while some excep- 
tionally luminous events can show their signature in polar- 
ized light. This could point out to the transient occurrence 
of QPOs which may appear repeatedly during the bright 
flares, which seems to be a rather natural possibility. This 
will be discussed in the next sections in more detail. 

3. Modeling 

In this section, we first describe in detail our emission 
model, which is mainly based on synchrotron emission from 
accelerated electrons in the inner parts of a relativistic ac- 



cretion disk. We also describe the ray-tracing method used 
and the predictions of the model. 

3.1. Fluctuations of the inner parts of an accretion disk: A 
possible description for the observed signal? 

3.1.1. Emission model 

Since the discovery of X-ray and NIR flares from Sgr A*, 
several flaring regions theories tried to describe the physics 
behind them, varying from abrupt changes in the accretion 
rate of a Keplerian disk (Melia et al. 2001) to the interaction 
of the accretion disk with nearby stars (Nayakshin et al. 
2004). Although none of these scenarios can be ruled out, 
there are observational evidences that give more support to 
some of them. 

For example, as we mentioned before, the frequently ob- 
served rate of NIR flares (four to flve flares/day) makes it 
hard for disk-star interaction or tidal capture scenarios to 
describe the events. Of special interest to us are the ob- 
served quasi-periodic flux modulation during the NIR and 
X-ray flares. The recent unambiguous discovery of (quasi- 
) periodicity in an active galaxy (RE J 1034+396) reported 
by Gierlihski et al. (2008) brings more support to the idea 
that the similarity in the behavior of black holes extends 
from stellar-mass black holes to super-massive ones. The 
most interesting scenario could be a relation to the orbital 
time scale of the accretion disk, with a possible connection 
to the plunging region which feeds the black hole through 
a channel of inflow or a possible clumpy infalling flow. As 
we describe here and in the next section, our interpreta- 
tion of these variable signals (which relates flux modula- 
tions mainly to lensing and boosting effects) can open a 
new window to study physics in very strong gravitational 
regimes, very close to the event horizon of black holes. 

NIR spectroscopy has shown that a power-law flt, ex 
(where Fj^, u and a are the flux, frequency and spectral 
index respectively), can describe the observed spectrum of 
NIR flares. Although all the observations agree with the fact 




Fig. 14. The geometry we considered in our emission 
model. The accretion disk around the black hole lies on 
the y — r plane and h is the unit normal vector of the disk. 
The direction of the magnetic fleld lines (B) correspond- 
ing to the disk frame is deflned by two angles r] and i/j. n 
represents the momentum of the emitted photon and its 
direction determined by a set of angles [(5e,0e]- The distant 
observer is looking into the system along a line of sight 
inclined by a certain angle i. 
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Fig. 15. Snapshots of orbiting anomalies inside the accretion disk as they appear to a distant observer looking along 
a line of sight inclined by 0° (relative to the normal to the disk), at times after |T, |T, |T and |T (left to right). 
Each row shows how the event evolves in time for different values of the characteristic shearing time scale: Tgh = oo 
(a-d), Tsh = 2.0 (e-h) and Tgh = 1-0 (i-1). The spin of the black hole is set to 0.5. The dotted and dashed lines indicate 
the position of the event horizon and marginally stable orbit, respectively. Each row has been scaled by its maximum 
intensity for illustrative purposes, a and (3 are the projections of the impact parameter of the emitted photons as the 
coordinates on the sky of the observer. Both coordinates are labeled in r^(— 5/ias) units. 



that NIR flares show a soft spectrum (a > 0), the value of 
the spectral index is still not well determined. The first NIR 
spectroscopy observations in July 2004 (Eckart et al. 2004) 
proposed the a value to be ~ 0.8 — 1.3 during the peak of 
the flare. In 2006, Gillessen et al. (2006) observed a correla- 
tion between flux and spectral index in their observations. 
However, recent observations by Hornstein et al. (2007) are 
consistent with a constant spectral index, a = 0.6± 0.2. 

The actual value of the spectral index shows its impor- 
tance in the modeling of the physical process responsible 
for the flaring emission. Some current models (Melia et al. 
2001, Liu et al. 2006, Yuan et al. 2007), predict that dur- 
ing flares a fraction of electrons near the event horizon of 
the black hole are accelerated. This can be described in the 
simplest form by a power law distribution in the energy of 
radiating electrons, A^(7) = A^oT"-^ where N{j)^ Nq^j and 
p are the energy distribution function of electrons, normal- 
izing constant, Lorentz factor of the electrons and the en- 
ergy spectral index respectively. For high values of a one 
will need a sharp cut-off to the energy spectrum of elec- 



trons (7c), while a lower value of a {a ^ 0.6) allows for a 
relatively milder distribution in the energy of electrons. Liu 
et al. (2006) have shown that simultaneous NIR and X-ray 
spectral measurements can constrain the parameters of the 
emitting region well. 

Before describing the details of our simulations, here 
we discuss how the existing observations limit the possi- 
ble range of free parameters. Observationally it is proven 
that in the Sgr A* spectrum a turn-over frequency in the 
sub- millimeter to NIR range exists. By using the turn-over 
frequency relation, Uc = 2.8 x 10^57^ GHz, one can put an 
upper limit on 7^5, where B is the magnetic field strength 
in Gauss. Here we have used 7c = 100 and B = 60G which 
give the best fit to the NIR/X-ray models that already exist 
(Liu et al. 2006; Eckart et al. 2008a). 

In our simulations we first considered a scenario in 
which the main flare is caused by a local perturbation of 
intensity close to the marginally stable orbit (via magnetic 
reconnect ions, stochastic acceleration of electrons due to 
MHD waves, magneto rotational instabilities (MRI) etc.). 
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Fig. 16. Snapshots of an orbiting anomaly inside the accretion disk as it appears to a distant observer looking along a 
line of sight inclined by 0° (a-d), 30° (e-h) and 60° (i-1) (relative to the normal to the disk) at times after \T, f T, 
and |T (left to right). See also the caption of Fig. 15 



These instabilities spread out and produce a temporary 
bright torus around the black hole. In this scenario, the 
mentioned variabilities are mainly due to relativistic flux 
modulations caused by the presence of an azimuthal asym- 
metry in the torus. 

Simulations are dealing with two important velocities: 
radial {vr) and azimuthal ('i;^). The radial velocity can be 

parameterized in the following way: Vj. ~ (^^f^)(^^)^ 
which depends on the ratio of the stress to the magnetic 
field energy density, jSj,^ and the ratio of the magnetic en- 
ergy density to the thermal pressure, pp (Melia 2007). The 
use of the typical values oi jSp and 13 from MHD simula- 
tions give us an estimation {PpPu ^ 0.05). This leads to a 
radial velocity of the order of 0.1 (^^). For the azimuthal 
velocity, we assumed that above the marginally stable or- 
bit the plasma is in a Keplerian orbit, v^k = ( ^ Z^^^'^^^ ) 

where A = — 2r — a^, and inside the plunging region the 
matter experiences free fall with the same angular momen- 
tum as at the marginally stable orbit. 

Furthermore, two important time scales are at work: 
heating and cooling time scales. The heating time scale 
strongly depends on the physical processes which act as 
the engine of the whole event (MHD instabilities, magnetic 



reconnect ions etc.), and one can just put an observational 
constraint on that according to the averaged observed rise 
time of the events {trise ^ 40min). Cooling time is mainly 
controlled by the Keplerian shearing and synchrotron loss 
time, isyn = 5 x 10^B~2i/-2 min, where B must be set in 
Gauss and u is in GHz. 

3.2. Hot spot model and fluctuations of the inner parts of 
the accretion disk 

Since all these physical processes happen very close to the 
black hole and in a very strong gravitational regime, we 
must take into account the effects of curved space-time. To 
simulate the changes in paths and polarization properties 
of photons from the emitting electrons to the observer at 
infinity, we have used the KY ray-tracing code (Karas et al. 
1992; Dovciak et al. 2004). KY is able to calculate ah the 
effects of GR, like light bending and changes in the emission 
angle, changes in the polarization angle of photons, gravita- 
tional lensing and redshift, Doppler boosting (since matter 
inside the accretion disk is in orbit) and frame dragging (in 
case of Kerr black holes) in a thin disk approximation. In 
the geometrical optics approximation, photons follow null 
geodesies, and their propagation is not affected by spin- 
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Fig. 17. Flux modulations of an evolving perturbation close 
to the marginally stable orbit of a Kerr black hole with spin 
parameter of 0.5. The light curves show how different values 
of shearing parameter and inclination affects the resultant 
light curve. 




Time m 

Fig. 18. Flux modulation (top), changes in polarization an- 
gle (middle) and degree (bottom) for a spot on a circular 
orbit at 1.0 x r^s (solid), 1.2 x rms (dotted), 1.5 x r^s 
(dashed), 2.0 X Tms (dott-dashed) around a Kerr black hole 
with a spin parameter of 0.5. The time unit is the orbital 
time scale (T). 



spin interaction with a rotating BH (Mashoon 1973). This 
means that wave fronts do not depend on the photon po- 
larization, and so the ray tracing through the curved space- 
time is adequate to determine the observed signals. Since 
our analysis is focused on high frequency regimes, we have 
mainly ignored radiative transfer effects. 




Fig. 19. Cross correlation between the changes in flux and 
polarization angle (top), and flux and polarization degree 
(bottom) for spot circular orbits at 1.0 X TjYig (solid), 1.2 X 
'^ms (dotted), 1.5 X Tfns (dashed), 2.0 x Vfjis 

(dott-dashed) 

around a Kerr black hole with a spin parameter of 0.5 (Fig. 
[is]). The time unit is the orbital time scale (T). 



To make KY work, we must initialize the properties of 
radiated photons at each point of the emitting region. The 
straightforward way is to define the intrinsic emission of 
each point in the context of the Stokes parameters. For the 
flux densities (mJy) and source sizes (/ias) of Sgr A*, op- 
tically thick synchrotron emission in the NIR can safely be 
disregarded (see the discussion in Eckart et al. 2009). Since 
in this report we focus only on the modeling of the NIR 
flares, it is sufficient to pick up a model for the energy dis- 
tribution of non-thermal electrons, radiating in an optically 
thin regime: 







7 > 7c 



which leads to the formulae for polarized emission: 



oc n (Bsin6le)^''^' V"'' 
= Ul cos (2Xe) I,^ 
= Ul sin (2xe) lu 



(7) 



(8) 
(9) 
(10) 

(11) 



where n, Xe, and He are the number density of the elec- 
trons, the angle between a reference direction and the plane 
of an observer co-moving with the disk frame and degree of 
linear and circular polarization, respectively. Throughout 
this paper we assume that in the NIR regime the light is 
not circularly polarized (lie = 0). lu^Qy^Uu, and rep- 
resent the Stokes parameters. Oe is the angle between the 
direction of the magnetic field and the direction toward the 
co-moving observer 



3(7?, ^e, ^e) = arccos( J -— ^ 



(12) 



where we have picked the disk co-moving frame as the ref- 
erence, so that the normal to the disk, n, coincides with 
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Fig. 20. A sketch showing how the pattern recognition coefficient (S) is defined. S is the multiphcation product of 
the cross correlation between observed and theoretical flux (cr/), observed and theoretical polarization angle (cr^) and 
observed and theoretical polarization degree (an) hght curves following Eq.[l6] Note that S(r') (g) is defined as a function 
of = r + (where r is the time lag in units of minutes) in order to make it easier to match the position of its peaks 
with the position of the lensing events in the observations. 



the z direction. Geometrical orientation of the dominant 
component of the global magnetic field vector (S) in this 
frame can be defined with a set of angles [t^, ?/^] (0 < 7^ < tt, 
< < 27r, see Fig. 14). (5e is the angle between the direc- 



tion of the photon momentum (/^) and the normal to the 
disk (n): 



= arccos ( 



(13) 



where v is the four-velocity of matter in the disc. 0e is the 
azimuthal emission angle which is defined as the angle be- 
tween the projection of the three momentum of the emitted 



photon into the equatorial plane and the radial tetrad vec- 
tor: 



(14) 



Having all the needed information at hand, by fixing 
free parameters like the spatial density distribution of the 
emitting electrons, magnetic field strength, flux spectral in- 
dex, and the global configuration of the magnetic field, one 
can simulate the expected images that a distant observer 
will measure for different possible inclinations with respect 
to the black hole/accretion disk system. Here we will show 
how compact azimuthal anomalies inside a uniform den- 
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Fig. 21. Pattern recognition coefficients as a function of time for our sample of observations [13 June 2004 (a), 30 July 
2005 (b), 1 June 2006 (c), 15 May 2007 (d), 17 May 2007 (e) and 28 May 2008 (f)]. Dotted lines show the median (red), 
68.3% (green), 99.8% (blue) and 99.99% (violet) significance levels derived from 10^ random red noise light curves. 



sity distribution of the emitting plasma can reproduce the 
observed behavior of Sgr A* in high frequency regimes. 

The synchrotron cooling time at NIR frequencies is only 
of the order of a few minutes, which is significantly shorter 
than the time scale of the observed flares. One way that a 
spot can survive long enough to be responsible for the ob- 
served ~ 4 — 5 time flux modulations is that a SSC mecha- 
nism up-scatters the sub-mm seed photons to the NIR and 
X-ray frequencies (Eckart et al. 2006a-c, 2008a). The other 
possibility is that the heating time of the orbiting NIR com- 
ponent is related to the rise time of the main flare, r oc trise 
where the emissivity profile of the emitting component fol- 
lows F{t; to, r) = Fo exp {-{t - tof/2T^). 

The above discussion demonstrates that it is essential 
to consider the gravitational shearing time scale as a vari- 
able in the simulations. We have implemented this effect in 
our modeling by introducing a dimensionless characteristic 
shearing time scale: 



Tsh = 



T{rspot) 



T{rspot + ^o) - T{rspot - ^o) 



(15) 



where for the initial spatial distribution of the relativistic 
electrons we have used a spherical Gaussian distribution 
with its maximum being located at the radius r = rspot 
with FWHM of ro = Ir^. In our simulations, Tsh varies 
between oo and 0.8. Tsh = oo corresponds to the situation 
in which the spot preserves its shape for a long time. The 
mechanism that stabilizes the spot is not known, although 
several possibilities have been proposed and explored in the 
literature. In particular, the existence of persistent vortices 
on accretion disks has been explored (Abramowicz et al. 



1992; Adams & Watkins 1995). However, our modeling sug- 
gests that shearing effects are indeed important and can 
be well represented within the multi-component scheme or 
in the spiral pattern scheme generalizing the original spot 
scenario. The Tgh = 0.8 corresponds to a pure Keplerian 
shearing for a spot located at the marginally stable orbit of 
an extreme spinning black hole (a = 1). One must note 
that MHD simulations are unable to produce hot spots 
with long life-times comparable to the observed flare time 
scales (Hawley et al. 2001, Krolik et al. 2002). Schnittman 
et al. (2006) used a model for the creation and annihila- 
tion of spots with short life times, distributed by random 
phase within a belt close to the marginally stable orbit. 
They show that their model can resemble the observed 
quasi-periodicity in the X-ray light curves. Eckart et al. 
(2008a) successfully modeled the simultaneous NIR/X-ray 
flares of Sgr A* by following the same basic idea of a multi- 
component model. As mentioned in §2, following the results 
of Pechacek et al. (2008) these kinds of multi-component 
models are able to reproduce the red noise behavior of PSDs 
while some bright individuals can show their signatures in 
polarized light (see Fig. 25 and also discussion in Eckart et 
al. 2008a). 

3.3. Results of the modeling 



Figure 15 shows how a hot spot is created and evolves in 
time for three different values of Tsh- A comparison between 
the rows shows how pure Keplerian shearing disrupts the 
initial shape of the spot and produces an elongated spiral 
shape. Figure [16] shows the apparent images of a spot with 
a mild shearing environment {tsh = 2.0) for three different 
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Fig. 22. A selected 100 minutes window of the white noise 
added simulated light curves of an orbiting spot plus a tem- 
porary variability in the accretion disk. The flux density is 
presented in arbitrary units (top). Middle and bottom pan- 
els show the changes in the angle (degree) of polarization. 
The spin of the black hole is set to be a = 0.5 and the 
inclination fixed on i — 60°. 



inclination angles. When we look face-on at the event {i o:^ 
0°), there are no modulations by the relativistic effects. For 
higher inclination angles {i = 30° and i = 60°) lensing 
and boosting effects play major roles in the observed flux 
modulations. Specially for high inclinations one can see how 
an Einstein arc develops when the spot passes behind the 
black hole and how photons coming from the accretion disk 
are blue-shifted on the left hand side of the image according 
to Doppler boosting. This set-up allows us to simulate light 
curves for a wide range of possible free parameters, mainly 
by covering the range of all possible inclinations and spins 
of the black hole. Figure [T7| shows examples of light curves 
for different values of inclination and shearing time scale. 



3.3.1. Pattern recognition analysis: signatures of lensing 
effects 

Figure [18] shows the typical magnification of flux, polariza- 
tion angle and degree of an orbiting spot emission as a func- 
tion of time. Here we showed the spots located at different 
distances from the black hole. These plots indicate the typ- 
ical behavior of light curves when the strong gravitational 
regime is prominent and the strong lensing and boosting is 
happening. As one can see, the sweep in the polarization an- 
gle precedes the peak in flux magnification, while the peak 
of the degree of polarization follows the magnification peak 
(see Broderick and Loeb 2006 for a detailed discussion). 




Fig. 23. Same as Figs.[3]and|5]for the simulated light curves 
of Fig. [22] The top shows the autocorrelation of the flux 
while the bottom shows the Lomb-Scargle periodogram. 
The colored regions indicate the position of the peaks cor- 
responding to the 0.8 — 1.2 X rmso orbital time scales of a 
Kerr black hole with spin 0.5. Dotted lines show the median 
(red), 68.3% (green) and 99.8% (blue) confidence levels of 
the red noise. 




Figure 19 shows more clearly this typical behavior accord- 
ing to the position of the corresponding cross-correlation's 



Fig. 24. Top (bottom): Cross correlation between the to- 
tal flux and polarization angle (degree) of the simulated 
light curve (Fig. [22|. Dotted lines indicate the correspond- 
ing cross correlation functions for the 30 July 2005 flare 
as the case most similar to the predictions of the hot spot 
model. 



peaks. It is particularly apparent that even with changing 
the position of the emitting source with respect to the black 
hole this effect remains more or less the same. 

This constant behavior encouraged us to check whether 
or not this type of pattern is manifested in our NIR light 
curves of Sgr A*. For this purpose we used a simple pattern 
recognition algorithm mainly via a multiplication of differ- 
ent cross correlation functions. Similar pattern recognition 



18 



Zamaninasab et al.: NIR flares of Sgr A* 



algorithms are used to identify gravitational wave signals 
from noisy data (Pappa et al. 2003, Goggin 2008). 

We derived a final pattern recognition coefficient prod- 
uct by multiplying two cross correlation functions (namely 
ai and (jj) as below: 



y^(<7i + l)(<T, + l) 



1 



(16) 



By applying this procedure for the cross correlation func- 
tions between the observed flux and the theoretical magnifi- 
cation light curve (cr/), cross correlation functions between 
observed polarization angle and the theoretical polarization 
angle light curve (cr^) and the same function for the degree 
of polarization (an), we derived a final pattern recognition 
coefficient product (H). Figure 20 shows how this pattern 



recognition coefficient is definecL 

Figure 21 shows the result of our pattern recognition 
analysis for the polarized flare events discussed in this pa- 
per. In order to estimate how significant the peaks in the S 
function are, we have repeated the same analysis for 10^ 
random red noise light curves simulated with the same 
method mentioned in Sect. 2.2.2. In all but one case the 



patterns shown in Fig. 21 can be identified at the > 3a to 
5(7 level. This shows that strong lensing patterns are signif- 
icantly manifested in our sample of NIR light curves. The 
strong lensing pattern detected is a further indicator for the 
existence of a (clumpy) accretion disk around Sgr A*. 

3.3.2. A spotted accretion disk? 



Figure 22 shows a selected 100 minutes window of the resul- 
tant light curves of the flux density, polarization degree and 
angle for a spot with constant shape {tsh = cxo), orbiting 
around a Kerr black hole (a = 0.5) close to its marginally 
stable orbit (r = 1.1 X TjYiso)' The line of sight is inclined by 
60 degrees {i = 60°). Gaussian white noise has been added 
to the simulated data in order to make the comparison of 
the periodicity and cross correlation with corresponding ob- 
servational results easier. The level of the noise and the av- 
erage error bars have been set from the average rms of the 
corresponding observed light curves. The surface brightness 
of the components follows a profile similar to F(t; to, r) with 
r = 25. The maximum degree of polarization which can be 
achieved via synchrotron mechanism is around 70%. Since 
any kind of deviation from the ideal isotropic distribution 
of electrons around the magnetic field lines will suppress 
the degree of linear polarization, we set the initial value 
for the radiation from the spot to be 50%. We assumed 
that the photons originating from the non-flaring part of 
the accretion disk are weakly polarized (~ 1%), since the 
main population of its NIR photons have thermal origin 
and relativistic electrons are randomly distributed around 
the magnetic field lines. 



Figure [23] shows the result from autocorrelation and 
Lomb-Scargle analysis (similar to Figs. |3] & |5|, and Fig. 
24 shows the results of the cross correlation analysis (sim- 
ilar to Fig. |8|. The peaks close to the minute time-lag 
are of specialinterest since they have the least dependency 
on the choice of free parameters and are mainly related to 
the basic idea that the flux modulations are caused by rel- 
ativistic effects. Figures 23 and 24 show that an orbiting 
spot is able to produce the same cross correlation pattern 
observed in our sample, but the orbital frequency of the 




Fig. 26. Same as Fig.[22jbut a selected 200 minutes-window 
of the multi component scenario. 




Fig. 27. Same as Figs. 23 for the multi-component scenario 

(Fig.m. 



spot will be detected to be mainly a quasi-periodical sig- 
nal. The corresponding correlation functions of the 30 July 
2005 observation are over-plotted in Fig. [24] for a better 
comparison. The main peak close to r = coincides very 
well for both observation and simulated cross-correlations. 

Furthermore, we have simulated a spotted accretion 
disk. In this case spots are born, evolve and finally fade 
away as a function of time. These anomalies are distributed 
in the inner part of the accretion disk in a belt between 



1 



Radial and azimuthal distribution of the 



anomalies are completely random and their distribution in 
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Fig. 25. Same as Fig. 16 for the multi-component scenario. The observer's inclination is set to 0° (a), 30° (b), 45° (c) 



and 70° (d) (see the discussion in Sect. 3.1 and Table. A. 2 in Eckart et al. 2008a). 
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Fig. 29. Snapshots of the orbiting spiral pattern as they appear to a distant observer looking along a line of sight inclined 



by 0° (a), 30° (b), 45° (c) and 70° (d) [relative to the normal to the disk]. See also the caption of Fig. 15 





servational results easier. As can be seen in Fig. 27 the 
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Fig. 

(Fig. 



28. Same as Figs. 24 for the multi-component scenario 
261). 



time follows a simple Poisson point process (see Pechachek 
et al. 2008). 

Figure [25] shows the snapshots of the spotted disk sce- 
nario as viewed by a distant observer from different incli- 
nation angles. The resulting magnification and polar imetric 
light curves are depicted in Fig. [26] White noise is added 
to all light curves in order to make the comparison between 
ZDCF and Lomb-Scargle results with the corresponding ob- 



random distribution of the spots strongly suppresses the 
periodic signal in the Lomb-Scargle periodogram, while the 
cross correlation between the magnification and the changes 
in the polarized flux still carries a significant signal from 
the modulations influenced by strong gravity (Fig. 28). As 
a conclusion one can say that a general relativistic simula- 
tion of turbulences in the inner parts of an accretion disk 
resembles the observed behavior of Sgr A* very well. 

One must note that for the simulated data, we have 
chosen typical values for spin and inclination (consistent 
with previous results by Eckart et al. 2006b and Meyer 
et al. 2006a,b), not trying to fit the actual data. As we 
can see, the data analysis and simulations show that even 
if there does exist a low-level activity physical process in 
Sgr A* which can be explained as random red-noise, the 
kind of observed correlations between the flux and polar- 
ization data cannot be produced via a completely random 
process (without taking into account relativistic effects). 
Actually, the similarity between Figs. 24 and [28] and what 
is already observed in at least three different Sgr A* flares 
(Figs. |5] and [8|, supports the idea of orbiting matter around 
the galactic super-massive black hole. 



3 A. Alternative models 

As mentioned before, the strength of ZDCF's peaks in Fig. 
[3] seems to be correlated with the brightness of flare events. 
This dependency on the brightness of the flares and also 
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the shape of the resulting ZDCFs are very similar to the ex- 
pected autocorrelation if the flux variations originate from 
geometric "light echoes" (Bursa et al. 2007; Fukumura et 
al. 2009). The expected periodicities from these models are 
different from the ones observed here, but their very close 
similarity keeps this issue open to future discussions. If a 
more detailed modeling of random flares happening at a 
distance from the black hole can show a possible ~ 20 min- 
utes QPO, then the light curve observed on 15 and 17 May 
2007 could be the first observation of second images created 
from separated photon "bunches" (Fukumura 2009). 

Furthermore, shocks in relativistic jets can produce 
a correlated total flux and polarized intensity. There is 
a possibility that under special circumstances (e.g spe- 
cial inclination, special magnetic field configuration etc.) a 
(episodic) relativistic outflow could produce the same corre- 
lation and pattern we have discovered in our observations. 



4. Geometry of the emitting region 

In this section we discuss the basic assumption of the ex- 
istence of an azimuthal asymmetry in the accretion flow 
of Sgr A*. Recently, Falanga et al. (2007) and Karas et al. 
(2007) discussed that a global spiral pattern of disturbance, 
with an orbiting speed not directly associated with the un- 
derlying Keplerian velocity, can fit the observed NIR and 
X-ray modulations of Sgr A*. Their model has been used 
to fit the observed X-ray flare on 31 August 2004 (Falanga 
et al. 2007). Here we will show how a combination of po- 
larimetric observations and the next generation of VLTI 
measurements can reveal the geometry of the emitting re- 
gion of these high frequency photons. 

In order to reproduce the same density profile in the 
inner part of the disk, we have used a spiral pattern, char- 




Fig. 30. Light curves of the normalized flux (top, black) 
and polarization angle (bottom, blue) observed at infinity, 
obtained form general relativistic modulations associated 
with a spiral pattern (left) or a compact spot (right). Light 
curves are shown for three different inclinations: i = 30° 
(solid line), i = 60° (dashed hue) and i = 80° (dotted 
line). The modulations happen twice as often for the spiral 
pattern than for the spot in the same interval, according 
to the existence of two symmetrical arms. Flux values are 
normalized according to their average values. The spin of 
the black hole is set to be zero (a = 0). 



acterized by the emissivity function given below: 
/j, (X sm^((j) + alog( — )) 

^0 



(17) 



while the power-law index, 7, describes the overall radial 
decrease in the emissivity, arctan(a) is the pitch angle and 
To determines the outer radius where the pattern fades 
away. This parametrization approximates a spiral pattern 
evolving in the background of a Keplerian disk. The main 
emission mechanism is assumed to be synchrotron radiation 
from relativistic electrons with the same energy distribution 
as Eq. (1). The shape of the spiral pattern highly depends 
on 7, /3, a and tq. We parameterize our model in such a way 
that it facilitates a straightforward comparison with Tagger 
et al. (2006). 7 = l,a = /3 = 6 and rp = 20 give the best 



approximation (see Fig. 29). Figure 30 shows the changes in 



the flux and polarimetric quantities, as measured by a dis- 
tant observer for different inclinations. As mentioned before 
(and discussed in detail by Falanga et al. 2007 and Karas 
et al. 2008), such a spiral pattern can produce the same 
typical behavior in flux modulation as that caused by an 
azimuthal anomaly in the accretion disk. 

The main idea in our simulations is to include the po- 
larimetric radiation transfer in curved space- time, which 
provides the possibility to compare the behavior of different 
geometrical configurations. It seems that even though both 
of these geometrical set-ups (spiral shape or compact spot) 
show the same behavior in flux, the observed polarized flux 
will behave significantly different. In order to reduce the ef- 
fects of our specific assumptions about the emission process 
and magnify the signatures which result only from differ- 
ent geometrical structures, we have chosen a toy model for 
the initial polarized emission: = I^^ JJu = = 0. 
Figure |30] shows that the spiral pattern is unable to pro- 
duce strong changes in the polarization angle, while a com- 
pact azimuthal source produces a highly variable polariza- 
tion angle. Due to abberation, photons that come from dif- 
ferent parts of the disk are polarized differently, even for 
Schwarzschild black holes. For the Kerr case, the rotation 
of the polarization vector will be added because of gravita- 
tional frame dragging. The dependency of these changes on 
the position of the emitted photon is depicted in Fig. [31 
The top panel of Fig. [31] shows how the polarization vector 
of the emitted photons will be rotated due to the strong 
gravity of the central black hole. As one can see, there is 
a clear knot visible in this contour graph. If the radiating 
source passes through this knot, the observed polarization 
angle will swing dramatically (for a detailed discussion see 
Dovciak et al. 2008). As the bottom panel of Fig. 31 shows, 
when a compact azimuthal source orbiting around the black 
hole is close to its marginally stable orbit, this nod will be 
passed. The amount of change in the polarization angle 
depends on the compactness of the source, its position rel- 
ative to the black hole, the inclination of the observer and 
the spin of the black hole. On the other hand, if the flare- 
emitting region is a deformed pattern extending in both 
radial and azimuthal directions, the swing in the polariza- 
tion angle will not be strong. It is because of the fact that 
in each point of the light curve the polarization angle is the 
average of photons that come from different parts of the 
disk and have a different polarization angle (Fig. [31]). As a 
result, the swings in the angles cancel each other out, and 
no significant swing will be observed. 
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Fig. 31. Top: Contour graphs showing the changes in 
the polarization angle due to general relativist ic effects. 
Bottom: Spatial emissivity distribution of a compact hot 
spot. The observer's inclination is z = 45° and the black 
hole is assumed to be spin-less. The observer is located on 
the top of the pictures. The innermost stable orbit is shown 
in both images (solid line in top image and dashed line in 
the bottom). The graphs are represented in the coordinates 
X = r cos (p^y = r sin cf) in the equatorial plane where r and 
(j) are Boyer-Lindquist coordinates. Units are Tg in the x-y 
coordinates. 



The presence of changes in the degree of the NIR polar- 
ization and the swings in the observed polarization angle 
simultaneous to the flux magnifications supports the idea 
that the geometrical shape of the sources is dominated by 
compact azimuthal asymmetries rather than by radially ex- 
tended spiral patterns (see also discussion in Sect. 6 about 
the differences in the centroid motions of these two types of 
geometries). Here we must note that as Tagger et al. (2006) 
have mentioned, a second initial configuration was used in 
their simulations (simulation no. 2, Tagger et al. 2006). In 
that set-up, a clump of matter starts spiraling towards the 
black hole and produces the spiral Rossby wave. There, the 
ratio of the surface brightness between the hot core mov- 
ing radially inward and the tail which is produced plays 
a critical role in the resulting light curves. If the hot core 
is dramatically brighter than the tail, there is no practi- 




Fig. 32. A sketch of the model (not to scale) of two extreme 
cases of global magnetic field configuration. Model I: the 
localized flare happens inside the accretion disk, where the 
dominant component of the magnetic fleld is toroidal {Bj). 
Model II: over-density of the accelerated electrons happens 
at the tip of a possible short jet (wind). Magnetic fleld 
lines {Bji) are assumed to be elongated toward the axial 
symmetry axis of the system, parallel to the normal to the 
disk (n.) 



cal difference between this scenario and the orbiting spot 
model. For a complete polarization study of Rossby wave 
instabilities, one needs simulated proflles of MHD surface 
densities as a function of time, which is beyond the scope 
of this paper. 



5. Magnetic field structure and geometrical 
orientation of the system 

In our simulations the structure of the magnetic fleld lines 
according to the black hole/accretion disk system can be 
controlled by a set of parameters [r^, ip]. Of special in- 
terest are two extreme conflgurations in which the enve- 
lope of the accelerated electrons is located inside the ac- 
cretion disk with a global toroidal magnetic fleld (model I: 
[77 = f 7 = f ]) oi" inside a region with magnetic fleld lines 
perpendicular to the disk (model II: [r] = O^ip = 0]). The 
latter can be interpreted as the spot being located at the 
tip of a jet with magnetic fleld lines aligned parallel to its 
symmetry axis (see Fig. 32). 



We assumed that the magnetic fleld inside the accretion 
disk is dominated by the toroidal component. This struc- 
ture is compatible with the results of several MHD simula- 
tions (Hawley & Balbus 1991; Hirose et al. 2004; De Villiers 
et al. 2003). However, in our case it is not clear if the cloud 
of accelerated electrons is located inside the accretion disk, 
somewhere above it in the corona, or even inside the tip 
of a possible short jet/wind. In fact, there is a possibility 
that a source structure in which an accretion disk is asso- 
ciated with a short jet/collimated outflow can explain the 
multi- wavelength behavior of Sgr A* (Markoff et al. 2001; 
Yuan et al. 2002; Eckart et al. 2005, 2006a, 2008a-c). Here 
we discuss how different orientations of the magnetic fleld 
lines inside the emitting region affect the resultant physical 
parameters that one can extract via fltting the model to 
the observed data. 
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X.red 1 






[deg] 


[deg] 




\j^mso\ 




[%] 


[deg] 


[deg] 




I 


0.5 


90 


90 


4 


1.1 


5 


50± 10 


60± 15 


15 ±20 


1.63 ± 1 


II 


0.5 








4 


1.1 


5 


50 ±10 


65=1=10 


110 ± 15 


2.95 ±1 


III 


0.5 


78.56 ± 12 


72± 18 


4 


1.1 


5 


60± 10 


50 ±20 


30 


1.54 ± 1 



Table 2. Final parameters resulting from the best Xred ^IR flare observed on 30 July 2005. 
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Fig. 34. The resulting orientations for the direction of pos- 
sible collimated outflow from Sgr A* (0) for model I (blue 
dashed-dotted line) and model II (green long dashed line), 
^-n^- [rninj ' Thc shadcd regions indicate the range of possible values 
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The KY code allows us to simulate light curves for a 
wide range of free parameters, so one can fit the model pa- 
rameters to the actual observed data (Meyer et al. 2006a, b, 
2007; Eckart et al. 2006a,b, 2008). For this purpose, simu- 
lations have been carried out, as they cover a wide enough 
range of possible inclinations (0.1° < i < 85°), the initial 
degree of linear polarization (0% < Ul ^ 70%), possible 
combinations of rj and ip, and different orientations of the 
whole system on the sky (0° < O < 90°). Here O defines 
the direction of the normal to the accretion disk projected 
on the sky. The spot has an initial radius of IRs and or- 
bits very closely to the marginally stable orbit of a Kerr 
black hole (a = 0.5). The dimensionless shearing time scale 
is fixed to be 2.0 {tsh = 2.0), and the ratio of the surface 



brightness of the spot to the torus is set to five. We have 
chosen these values according to existing results of several 
fits to the NIR flares (Meyer et al. 2006a,b; Eckart et al. 
2006b, 2008). 

As a first step, we focus on model I ([77 = ^^ip = ^]). 
In order to find the fit with best Xred h.diYe carried out 
a grid search in a (z — — Hl) parameter space. The steps 
for i, 6 and Ul have been chosen to be 5°, 5° and 5%, 
respectively. The least xled value was achieved for a high 
inclination angle {i = 60°), O = 15° and highly polarized 
source (H^ = 50%) (see Fig. [33] and Table 2). The grid 
search has been performed on 11^, i and 9 for the models 
I and II. The free parameters in model III are 11/,, z, 77 
and ip. Figure [34] shows the resultant possible range for the 
orientation of the correlated outflow/wind in this scenario. 
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Interestingly, the observed position of the mini-cavity lies 
well within this interval. Muzic et al. (2007) have shown 
how a collimated outflow from the position of Sgr A* can 
describe the observed motion of filamentary structures and 
the mini-cavity, according to their observations in the NIR 
L'-band. This observational evidence supports the idea that 
the NIR flares actually originate in accelerated electrons 
within an accretion disk. 

For model II ([/^ = 0, = 0]), the envelope of relativistic 
electrons has been located somewhere in the bottom of an 
outflow. The least Xred achieved for this model and related 
parameters are presented in Table 2. Possible directions of 
this jet structure are depicted in Fig. [34) As one can see. 



Xred values and observational facts support the idea that 
the NIR photons originate in the inner parts of an accretion 
disk rather than in an outflow (Thorne 1994, Blandford 
2001). 

One can also invert the question and try to constrain 
the magnetic field structure of the emitting region by fixing 
the direction of the outflow in the plane of the sky (model 
III). For this purpose, a Xred search has been carried out, 
assuming that the normal to the disk is fixed at O = 30°. 
This is the best value that can fit the observed position of 
the mini-cavity (Muzic et al. 2007). By a good approxima- 
tion, the best fit is achieved for a toroidal magnetic field 
structure (see Table 2). This shows that if any future ob- 
servations reveal traces of an outflow from Sgr A*, NIR 
polarimetry could be used as a tool to constrain its mag- 
netic field structure. 



6. Next generation of NIR interferometry 

According to its proximity, Sgr A* is known as the best 
candidate for studying the details of physical processes on 
the event horizon scales (Falcke et al. 2000; Bromley et al. 
2001). VLT and Keck telescopes have already achieved an- 
gular resolutions as high as 40mas in their imaging mode 
(Genzel et al. 2003; Eckart et al. 2004, 2006a,b; Ghez et 
al. 2005). In interferometric mode these telescopes can ob- 
serve with resolutions of the order of only a few mas. 
The next generation of VLT Interferometer instrument, 
namely GRAVITY, would be able to achieve resolutions 
of the order of 10/ias in its phase-reference tracking mode. 
GRAVITY will try to measure the position of the lumi- 
center of the flares with respect to a reference source. For 
such a measurement the la uncertainty is ±10/ias with an 
integration time of five minutes on an assumed K = 15 
point source. The real necessary integration time for that 
precision is of course directly dependent on the magnitude 
of the flare. This resolution is high enough to track the 
centroid motion of NIR images of Sgr A* with a resolution 
of the order of one Schwarzschild radius (Eisenhauer et al. 
2005, 2008; Gillessen et al. 2006). These unique features 
make the next generation of VLTI a perfect instrument for 
high accuracy astrometry of the matter around Sgr A*. 
Recently, Doeleman et al. (2008) proposed a VLBI configu- 
ration with a high enough time resolution to possibly track 
the centroid motion of sub-mm images of Sgr A* with even 
higher angular resolution. In this section we present the re- 
sults of our study on centroid motion of the NIR image of 
Sgr A*. 

Some predictions of the centroid motions related to the 
hot spot model have already been studied by several au- 
thors (Broderick & Loeb 2006a, b; Paumard et al. 2006; 



Zamaninasab et al. 2008b; Hamaus et al. 2009). In previous 
works it has been mainly assumed that the spots preserve 
their shape in time. This assumption has been based mostly 
on the presence of strong magnetic field lines or "unknown 
mechanisms" (Hamaus et al. 2009). We will show here that 
it is actually very crucial that the effects of gravitational 
shearing inside the accretion disk are taken into account, 
since the predicted results can change dramatically. The 
other very important fact to be considered in the simula- 
tions is due to the NIR photons which come from the body 
of the accretion disk and confuse with the spot's emission. 
Since GRAVITY can only achieve its 10/ias resolution in 
the phase-reference tracking mode and not in the imaging 
mode, even a constant confusion from the (thermal) elec- 
trons inside the accretion disk can not be reduced. 



Figure [35] represents how these two parameters men- 
tioned can affect the expected centroid tracks and may 
present a complication in detecting the plasma structure 
close to the event horizon of Sgr A*. The simulated paths 
in Fig. 35 belong to a configuration where the observer is 
aligned approximately face-on {i 0°). This is the incli- 
nation in which the clearest wobbling of the centroid is 
expected to be seen. Fig. 35 1 shows that the position wan- 
ders of the center of the images are large enough to be 
detected easily if GRAVITY could achieve its five minutes 
time-resolution (simulated points). Panels b & c in Fig. 35 
show how gravitational shearing changes the centroid tracET 
One can see in particular that in the case of pure Keplerian 
shearing the change in the center is detectable for approxi- 
mately only six points. The second row (Fig.|35j panels d-j) 
shows the result of the simulated tracks including the effect 
of confusion due to the NIR photons coming from the inner 
parts of the accretion disk. 

Figure [36] shows the centroid paths of spots orbiting in a 
mild shearing environment {tsh = 2.0) for three different in- 
clinations {i = 0.1°, 30°, 60°) in presence of confusion from 
the torus. Panel c shows that the detection of any position 
wander becomes even more difficult in the case of a shearing 
spot observed at high inclinations. The bottom row in Fig. 
[36l shows the centroid motion resultant from simulations 



including several emitting components (Fig. 25). The pho- 
tons originated from the spot are highly polarized, while 
the radiation from the torus is weakly polarized. This fact 
may lead us to a solution for compensating for this kind 
of confusion. Figure [37] shows how different the apparent 
images of the flares are when they are observed in total 
flux or two orthogonal polarized channels. One sees that 
the different centroid paths can be revealed if GRAVITY 
could achieve its 10/ias resolution in the polarimetry mode 
FigJU 

To summarize the results of this section, one can say 
that a non-detection of a position wander in NIR images 
of Sgr A* will not rule out a spotted disk scenario. On the 
other hand, a clear detection of such a position wander by 
GRAVITY would lead us to strongly favor the spot model 
and would open a new era in studying the physics of high 
gravitational regimes. 

7. Summary and conclusion 

We used a sample of NIR flares of Sgr A* observed in po- 
larimetry mode to study the nature of the observed vari- 
ability. Using the z-transformed discrete correlation func- 
tion algorithm, we found a significant correlation between 
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Fig. 35. The centroid motion of NIR images viewed from 0° above the orbital plane for a Kerr black hole with a spin of 
0.5. The first row shows the paths for the background-subtracted images for different values of the gravitational shearing 
time scale. Each column shows Tsh values oo, 2.0 and 1.0 from right to left. The second row shows the same paths for 
the images without the background subtraction. The error bars show the simulated positions as they are expected to 
be observed by GRAVITY (FWHM ~ 8.5/ias) with a time resolution of five minutes. The whole paths belong to three 
orbital times. A triangle indicates the beginning and a square the end of the track. Both coordinates are labeled in /ias 
units. The dotted and dashed lines indicate the position of the event horizon and the marginally stable orbit respectively. 



changes in the measured polarimetric data and total flux 
densities. This provides evidence that the variations prob- 
ably originate from inner parts of an accretion disk while a 
strong gravity's eflFects are manifested inside them. 

In order to obtain this information polarization data 
is indispensable since the corresponding signals are very 
difficult or even impossible to extract from simple total in- 
tensity light curves only (Do et al. 2008). The presence of 
significant signals from orbiting matter then calls for de- 
tailed modeling of the NIR polarized light curves in order 
to analyze the distribution of the emitting material and 
the magnetic field structure within the accretion disk (e.g. 
Eckart et al. 2006a, 2008a, Meyer et al. 2006ab,2007). 

In addition, the pattern recognition algorithm we em- 
ployed in this paper is an eflftcient tool to search for fiare 
events that carry the signature of strong gravity in light 
curves which are significantly longer than the orbiting time 
scale over which such an event can typically occur. Other 
methods that make simultaneous use of the entire light 
curve (like e.g. the Lomb-Scargle algorithm) tend to di- 
lute these signals and lower the possibility for significant 
detections dramatically (Do et al. 2008, Meyer et al. 2009). 



In order to constrain the physical properties of the emit- 
ting region we employed a relativistic disk model with az- 
imuthal over-densities of relativistic electrons. A combina- 
tion of a synchrotron mechanism and relativistic amplifi- 
cations allows us to fit the real observed data and make 
predictions about astrometric parameters of the accretion 
disk around Sgr A*. The modeled light curves show the 
same correlation between the fiux and polarimetric data as 
the one deduced from observations. 

Simulations have been carried out in a way that they 
cover a wide range of parameters, including the effects of 
gravitational shearing inside the accretion disk, the heating 
and cooling time scales, the inclination and the spin of the 
black hole. It is discussed how the observed swings in the 
polarization angle support the idea of a compact source 
for the emission, instead of radially extended spiral shapes. 
Furthermore, we present a model in which the observed 
NIR polarization angle can lead to confining the expected 
region for the expected outfiow/wind from Sgr A*. The 
model also predicts that when observations will be able to 
resolve the position of such an outfiow, the magnetic field 
structure inside the accretion disk could be confined. 
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Fig. 36. Similar to Fig. 
i = 30°, i = 60°). {di-cyi. 
mult i- component model (see Fig. 25). 



where different columns represent different values of inclination (from left to right: i = 0.1°, 
rows the centroid track of an orbiting spot with Tsh = 2.0. (d-f) shows the centroid for a 




Fig. 37. Apparent images of a flare event for a distant observer looking along a line of sight inclined by 45° in total flux 
(a), ordinary (b) and extra-ordinary polarized channels (c). 



Finally, the centroid paths of the NIR images are dis- 
cussed. In comparison with the results by Broderick & Loeb 
(2006a,b) and Hamaus et al. (2009), we have shown that 
the geometrical structure of the emitting region (elonga- 
tion of the spot according to the Keplerian shearing, multi- 
component structures, spiral arms, confusion caused by the 



radiation from the hot torus) can affect the expected cen- 
troid tracks. While all the mentioned geometries are able 
to fit the observed fluxes, we show how the future NIR 
interferometer GRAVITY on the VLT can break these de- 
generacies. The results of simulations propose that focusing 
GRAVITY observations on the polarimetry mode could re- 
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Fig. 38. Centroid paths for the same inchnations as Fig. 37 in the total flux (a), ordinary (b) and extra-ordinary (c) 
polarized channels. See also caption of Fig. 35 



veal a clear centriod track of the spot(s). We conclude that 
even though a non-detection of centroid shifts can not rule 
out the mult i- components model or spiral arms scenarios, 
a clear position wander in the center of NIR images during 
the flares will support the idea of bright long-lived spots or- 
biting occasionally around the central black hole. This pos- 
sible detection opens a new window for testing the physics 
very close to the edge of infinity. 
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